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Abstract We present a comprehensive review of the current state of fracture phe- 
nomena in transient networks, a wide class of viscoelastic fluids. We will first define 
what is a fracture in a complex fluid, and recall the main structural and rheological 
properties of transient networks. Secondly, we review experimental reports on frac- 
tures of transient networks in several configurations: shear-induced fractures, frac- 
tures in Hele-Shaw cells and fracture in extensional geometries (filament stretching 
rheometry and pendant drop experiments), including fracture propagation. The tenta- 
tive extension of the concepts of brittleness and ductility to the fracture mechanisms 
in transient networks is also discussed. Finally, the different and apparently contra- 
dictory theoretical approaches developed to interpret fracture nucleation will be ad- 
dressed and confronted to experimental results. Rationalized criteria to discriminate 
the relevance of these different models will be proposed. 
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1 Introduction 



The ability of viscoelastic fluids to fracture has been recognized in the sixties | Hut- 
|ton| ([T963 )] but remains much less documented that the breakdown of solid materials. 
The starting point for discussing fracture in viscoelastic fluids is coming up with a 
rigorous definition of fracture. The layman's definition of failure proposed for solid 



materials by Buehler |Buehler (2010 1] - failure occurs when the load bearing capac- 
ity of the material under consideration is significantly reduced or completely lost due 
to a sudden development - is not applicable for fluids, which accommodate arbitrar- 
ily large deformations after a finite time. At a more fundamental level Buehler argues 
that fracture of a material due to mechanical deformation can be understood as con- 
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version of elastic energy into breaking of chemical bonds or heat. This definition is 
more appropriate for fluids but does not allow to discriminate clearly the fragmenta- 
tion of liquids due to some hydrodynamic instability |Eggers and Vill ermaux] ((2008 )] 
and the fracture of fluids reminiscent of the fracture of solids. For instance, liquid jets 
serve as a paradigm for hydrodynamic instability leading to drop breakup through the 
Rayleigh-Plateau capillary instability: in this case, intermolecular (elastic) energy is 
converted at the surface vicinity alone. A tentative definition for the fracture of fluid 
could be the following:" the fracture of a fluid due to mechanical deformation, includ- 
ing flow, can be understood as a dissipation of bulk elastic energy into breaking of 
physical or chemical bonds." This definition excludes the hydrodynamic instabilities 
observed also in inviscid or Newtonian fluids of the field of fluid's fracture because 
elastic energy that come in play is of interfacial nature solely. 

Complex fluids denotes a (too) huge class of condensed-phase materials that 
posses mechanical properties intermediate between ordinary liquids and ordinary 
solids |Larson| ( [l999| l]. Many of them are "solids" at short time and "liquid "at long 
time, hence they are viscoelastic; the characteristic time required of them to change 
from "solid" to "liquid" varies from fractions of a second to hours. Numerous exten- 
sive studies on the rupture of entangled polymer melts in extension have been carried 
out in the past forty years. This allows to distinguish two zones of either viscoelas- 
tic rupture or elastic fracture in the Malkin-Petrie master curve intended to illustrate 
the different extensional responses with increasing strain rates of entangled polymer 
melts [Malkin and Petrie ( 1997|l]. These two types of rupture have been recently re- 



considered by Wang |Wang and Wang (2010)]: a yield-to-rupture failure transition 



is observed: the yield failure (called also ductile failure by other authors |Ide and 



|White| ( fT976| 1 1977| 1 1978| >] should be due to the yielding of entanglements, whereas 
at higher extensional strain rates, the specimen breaks up in a purely elastic regime 
(without any flow). 

Among visoelastic fluids, self-assembled transient networks, that consists of re- 
versibly cross-linked polymers in solutions constitute model systems for the physi- 
cist, with well defined structural properties and simple linear rheological behavior. 
The aim of this review is to embrace the current aspects of the phenomenology of 
transient network's fractures in a compressive manner. We first present what are tran- 
sient networks: systems and rheology. Second, we report on fracture's experiments 
and computer simulations of transient networks in several configurations: shear ge- 
ometry, Hele-Shaw cells and extensional geometry. This also includes, fracture prop- 
agation and the extension of the concepts of brittleness and ductility for viscoelastic 
fluids. Then, we describe the several theoretical approaches developed to describe 
crack nucleation in transient networks, and how they compare to experiments, before 
to conclude. 



2 Transient networks: systems and rheology 

Self-assembled transient (often called physical) networks are a class of complex ma- 
terials forming spontaneously 3D networks at thermodynamical equilibrium, that 
can transiently transmit elastic stresses over macroscopic distances. Transient self- 
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assembled networks are common in both natural and synthetic materials. They con- 
sist of reversibly cross-linked polymers in which weak interactions such as hydrogen 
bonds, hydrophobic interactions, van der Waals forces, or electrostatic interactions 
are responsible for cross-links formation. Because of the transient character of the 
junctions, and so of the thermodynamics equilibrium state of these systems at rest, 
they don't exhibit any aging nor yield stress contrary to other classes of soft out-of- 
thermodynamic equilibrium viscoelastic materials like dense particulate suspensions 
(colloidal glasses), or thermoreversible gels [Larson ( 1999)]. One of the major issue 
of transient polymer networks is to convey useful rheological properties to solutions, 
such as increased viscosity, gelation, shear-thinning or shear-thickening. They can 
be used as controlled drug delivery systems | Sutter et al (2007 1], rheological reg- 
ulators in polymer blends | |Kim et al] ( |2004 i], coatings, food, and cosmetics, or as 
matrix materials for tissue engineering [Kim and Mooney (1998)]. Self-assembled 
transient networks consist mostly of binary solutions of associative polymers, or 
ternary solutions of associative polymers and self-assembled surfactant aggregates 
even if supramolecular reversible networks formed by mixtures of small molecules 
associating by directional interactions has been also reported [Cordier et al (2008)]. 
Associating polymers are macromolecules with a part that is soluble in a selective 
solvent (often water), the so-called backbone or spacer to which two or more moi- 
eties that do not dissolve in this solvent, the stickers, are attached. The stickers may 
be randomly distributed along the backbone or may be grouped in blocks, the associa- 
tion of such polymers in solution has been studied extensively, and many reviews can 
be found in the literature, see for exa mple ||Larson|(|1999|); Winnik and Yekta| ( |1997| l; 
|Berret et al| ( |1997[ ); peng and Russel| ( |2006| ); |Chassenieux et al| ( |201 1} ]. Note that the 
restrictive definition of transient networks that we propose, excludes simple entan- 
gled normal or living polymer solutions and melts, where entanglements of polymer 



chains play the role of transient cross-links and exhibit elastic yielding |Malkin and 



|Petrie| ( fT997| ); |Boukany et al] ( |2009| )]. Telechelic polymers are often used as model 
linkers because they are architecturally simple: they consist of a long solvophilic 
mid-block with each end terminated by a solvophobic short chain (a sticker) JSe 



menov et al (1995)]. The stickers incorporate into the solvophobic domains of the 



aggregates and can bridge them via their solvent-soluble mid-block resulting in an 
attractive interaction between the aggregates. . The nature and the morphologies of 
the aggregates forming the network's junctions are versatile: (i) telechelic polymers 
in binary solution that self-assemble spontaneously into non-interacting flower-like 
micelles at low concentration and form three dimensional networks above a perco- 
lation concentration ]Annable etall ( |1993| ); [Serero et al| ( (2000| ); |Werten et al|p009] ); 

(2006)], (ii) surfactant vesicles 1 Lee et al ( 2005)1], (iii) lyotropic l amellar 



Seitz et al 



phases |Warriner et al ( 1997[)] (iv) wormli ke micelles | Ramos and Ligo ure| ( |2007) ); 



Lodge et al| ((20071); [N akaya-Yaegashi et al (2008); [Tixier et al| ( |20T0)i], (v) spherical 
micelles [Appell et al|(|1998[);|Tixier et al|(|2010])]7 (vi) oil-in-water [B agger Jorgensen 



et al ( 1997|);|Filali et al|(|1999[)] or (vii) water-in-oil microemulsion droplets [Oden- 



wald et al (1995)], (viii) photocrosslinkable nano-emulsions [Helgeson et al (2012)]. 



The rheology of transient networks is determined by the amount and the life time 
of the bridges. In the simplest case when the bridging chains are flexible and un- 
tangled, each bridging chain contributes about 1 kgT per unit volume to the elas- 
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tic modulus Go according to the theory of rubber elasticity [Green and Tobolsky] 
7946}; Tanaka and Edwards ( 1992ei; Yamamoto ( 1956| l]. The typical value of the 



shear modulus for these networks ranges between few Pa to several ten thousands 
Pa. The elastic response relaxes when the solvophobic blocks escape from the core. 
Often the escape is characterized by a single relaxation time so that the terminal re- 
laxation of the viscoelastic properties is characterized by a single Maxwell process 
well separated from the faster internal modes that characterize the conformational re- 
laxation of the chains | |Annable et aT| ( [T993] l; [Serero et al|pOOO] >; |Michel et al|p000] l; 



Filali et al (2001 1; Tabuteau et al|(|2009a)>; [Hough and Ou-Yang|(|2006[)]. Far above 



the percolation concentration, the viscoelastic relaxation time T is related to the av- 
erage lifetime of a connection | Green and Tobolsky (1946 )] which in turn depends 
on the breakage probability of a cross-link. The average lifetime of bridge will de- 
pend on its chemical nature, the external conditions and the physical state of the 



cross-links. It can vary through few ms |Tixier et al (2010 1], few seconds | Michel 



et~aT| (pOOO) ; |Filali et al| ([2001) ], up to hours ISerero et al| ( |2000| ; |Skrzesze wska 
et al| ( |2010p ; |Seitz et al (2006i]. The simplicity of the linear viscoelastic behavior 
of most of self-assembled transient networks is in contrast with their highly complex 
non-linear response that can vary widely | |Pe liens et al ( 2004b )] . Among all systems, 
the experimental model transient network consisting of oil-in-water microemulsion 



droplets reversibly linked by telechelic polymers | Filali et al (2001); Michel et al 



(2000)] is perhaps the unique one which exhibits the steady shear flow curve of a 
pure Maxwell fluid for both the shear stress and the first normal stress difference un- 
til it breaks | Tabuteau et al (2009a)]. In most other systems, the steady shear viscosity 
exhibits three flow regimes. Above a Newtonian plateau (linear regime), the viscosity 
increases considerably (shear thickening) prior to the onset of shear thinning at higher 
shear rates [Annable et a T|( fT9l^ ; |Xul^ ( [T9l^ ; |0^ and Serero| 



(|2001|);|P"eTiens et al|([2004a|);|Tripathi et al|(|2006[)]. However, the shear thickening re- 



gion is not always present [Michel et al (2000|);|Tixier et al|(|2010|);[Tirtaatmajda and| 



Jenkins (1997i; |Mewis et al| (|200r); Skrzeszewska et al (2010)]. A large number of 



constitutive theoretical models are based on the temporary-network kinetic model for 
telechelic polymers network theory [Tanaka and Edwards (1992e)] by applying the 
original ideas formulated by Green and Tobolsky [Green and Tobolsky (1946)] and 
Yamamoto [Yamamoto (1956)] and have been developed to tentatively capture the 



main features of non-linear rheological properties of these networks JMarrucci et al 
( 1993[); Ahn and Osaki ( 1995[); van den Brule and Hoogerbrugge (2000|); |Vaccaro and 
Marrucci|(|2000|);|Tripathi et al (|2006|);|Hernandez-Cifre et al|i|2003}L 



3 Fracture experiments 

Several geometries have been considered to explore fracture mechanisms in tran- 



sient networks: elongational flows using either extensional rheometer [Tripathi et al 



( 2006) )], pendant drop experiment s [Tabuteau et al| (|2009a| |201 ![)], shear fl ows in 
rheometric cells [Be rret and Serero| ( |2001| ); |Tabuteau et al| ( |2009a| ); |Tixier et al| ( 2010] ); 



Skrzeszewska et al] ( |2010 



Ignes-Mullol et al| ( |1995| ); 



i] or flows in Hele-Shaw cells [Zhao and Maher U993); 



Vlad et al ( 1999 ); Mora and Manha| ( 2010 1] . 
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Fig. 1 From Tabuteau et al 

(o N = <5;z 



2009a ]: (Left) shear stress (cy.) and first normal stresses difference 
) versus shear rate for a bridged microemulsion with a shear modulus Go = 1210 Pa and 



a relaxation X = 0.8 s. The continuous lines are fits (Maxwell model). (Right) Development of a fracture, 
occurring at the surface of the sample, for a shear rate equal to 0.9s~', corresponding to a critical first 
normal stresses difference Ojv = 1 190 Pa. 



3.1 Shear-induced fractures 



3.1.1 Stationary shear rate 



When submitted to a constant shear rate y, the measured shear stress ff™ of "brittle" 



transient networks first increases smoothly with the shear rate | Molino et al| ( 2000|; 
Berret and SerCTo| | |200T) ; |Tabuteau et al| ( |2009a| i; [Skrzeszewska et al|p01o| ; |Tixier| 



et al 



(2010 1] . However, above a critical shear rate y ~ z , where T is the relaxation 



time of the network, the flow curve exhibits a sharp discontinuity (Figure [T| and 
the viscosity decreases abruptly. Below this value, the branch of the flow curve can 
be Newtonian and shear thickening [ Berret and Serero] ( [20 01 1], or only Newtonian 
[Moling et al] ( |2QQQ> ; |Tabuteau et al f2009a| l; |Skrzeszewska et all ( [20TUl i], or New- 
tonian and shear thinning | Tixier et al ( 2010| l] depending on the experimental sys- 
tem under consideration. Note that for "brittle "transient networks, the first normal 
stresses difference Ni = a xx — oy Y drops suddenly for the same critical shear stress 
| Tabuteau et al (2009a)]. Authors of |Molino et al (2000)] were the first to suggest 
that the sudden drop of the stress in the flow curve of a transient network is the sign 
of a fracture propagation. 

However, the first unambiguous demonstration of shear-induced fractures in tran- 
sient networks was done in [Berret ancT Serero (2001 )], by using a flow visualization 
technique in a plate-plate transparent shearing cell. At low shear rate, the velocity 
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profile is homogeneous: the velocity decreases linearly from the rotating wall to the 
stationary one. Above the critical shear rate, the stationary velocity field within the 
gap of the cell exhibits a discontinuity, that defines a zone of fracture as the part 
of the fluid submitted to a high shear rate (~ 10 times the applied rate). This has 
been confirmed by Skrzeszewska et al | Skrzeszewska et al (2010)] by Particle Im- 
age Velocimetry: the fracture zone has an irregular shape and is rather wide on the 
order of a few hundred micrometers. With increasing overall shear rate, the width 
of the fracture zone increases. The fracture zone can happen everywhere in the gap 
and is different for every experiment. For non-adhesive gels, the fracture zone occurs 



generally at one of the wall |Berret and Serero (2001 1], and so appears as sliding, 
that is observed in other classes of complex fluids like pastes or concentrated emul- 
sions [Me eker et al| ( |2004 i], whereas for adhesive gels [Skrzeszew ska et al| ( |2010[ )], 
or rough wall surfaces |Tixier et al ( 2010) )], the fracture occurs in the bulk. A direct 
optical observation of the shear fractures | Tabuteau et al ( 2009a l] shows that above 



the critical stress, cracks open up all around the sample and grow rapidly. It is worth 
noting that the fractures are tilted 45° from the shear plane, perpendicular to the di- 
rection of the maximum extension (Figure [T}. Interestingly, for two very different 
systems |Tabut eau et al| ( f2"0"09a| ); |Skrzeszewska et al (2010)] that have perfectly New- 
tonian flow curve, before fracture occurs, the fracture shear stress in the steady state 
flow curve scales as (7„ ~ Go, where Go is the shear modulus of the network. In con- 
trast to solids, here the fractures heal over rapidly after the shear rate is switched off, 
and a new experiment can be performed after a few minutes with quantitatively the 
same behavior [Tabuteau et al (2009a); Skrzeszewska et al (2010)]. 



Above the critical fracture shear rate, strong fluctuations of the shear stress have 
been observed JSprakel et aT] ( |2009a|l;|Tixier et al|(p0T0) l; |Ramos"etlIl ( [20TT) ], that can 
lead to an apparent shear plateau flSprakel et al| ( |2009a] >] reminiscent of shear-banded 
flows observed in a wide variety of soft materials such as solutions of entangled 
wormlike micelles, colloidal suspensions and entangled polymer solutions | Fi elding| 
( 2007| ); |O"Im sted (2008)]. Such strong fluctuations of the shear stress have also been 
observed in solutions of entangled wormlike micelles at high shear rate and have 
been associated also with a rupture-like behavior as evidenced by Particle Tracking 
Velocimetry [Boukany and Wang (2007 1]. By using a novel class of transient net- 
works made of surfactant micelles of tunable morphology [Tixier et al (2010)] (from 
spheres to rods to flexible worms) linked by telechelic polymers, and coupling rhe- 
ology and time-resolved structural measurements, Ramos and Ligoure | Ramos et al 



( 201 1| >] clearly show that true shear-banding is not associated with strong fluctua- 
tions of the shear stress. Indeed, fluctuations of the shear stress are entirely correlated 
to the fluctuations of the degree of alignment of the micelles that can probe a frac- 
ture process. Sprakel et al [Sp rakel et al| ( |2009a| ] argue that the intermittent behavior 
observed in the stress response is due to repeated microfracture-healing events in 
the material. The cumulative distribution of the total stress drops A a during a frac- 
ture displays a characteristic power-law behavior, &(> Aa)°^ S<T~° ,8S , with & the 
stress probability distribution. The exponent is close to the value of 0.8 reported for 
true stick-slip motion [Feder and Feder| ( |1991) ]. Attempts to explain such scaling 
behavior involve the concept of self -organized criticality. 
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Fig. 2 From Sprakel et al ( 2009b ] : Transient shear stress response in a shear start-up run. Bottom panel 
shows the corresponding evolution of the velocity profile, in which each line represents the local fluid 
velocity (v x (y,t)) at a given position in the gradient direction. Snapshots of the simulation box (flow from 
left-right with a velocity gradient from top-bottom) illustrate the homogeneous initial configuration (left) 
and the final fractured state (right). 



Responsive particle dynamics simulations | Sprakel et al ( 2009b)] have also shown 
shear-induced fractures (Figure [2]i in transient polymer networks. Moreover these 
simulations have revealed a transition from shear banding to fracture upon increas- 
ing the overall polymer concentration, and emphasize the difficulty to define an un- 
ambiguous criterion to discriminate between banding and fracture. In the literature, 
banding is usually defined as situations where the velocity profile in the gradient di- 
rection is continuous but kinked, whereas it is discontinuous at the plane of fracture. 
It is worth mentioning than the the issue of velocity change within the gap is in fact 
more complicated [Manneville (2008 )] in addition to "classic banding" coexistence 
of yielded and unyielded regions with finite and zero rates/velocities have been ob- 
served in systems like colloidal glasses |Besseling et al ( 2010) !]. Unfortunately the 
spatial accuracy is typically limited to the micrometer range and cannot allow to dis- 
tinguish between narrow high shear bands and fracture planes. However the authors 
of | Sprakel et al (2009b)] show that a discontinuity in the variation of the first normal 
stresses difference N\ with the shear-rate is an unambiguous signature of a fracture, 
since normal forces should be largely reduced as all connections between the two 
shear bands across a fracture plane are broken. The discontinuity of N\ at the thresh- 



old has been observed experimentally for fracturing transient networks [Tabuteau 
et al ( 2009a| ); [Tixier et al] ( |2010| )], in contrast with shear-banding regimes |Tixier et al 
( 2010| )], where N\ increases with y as expected theoretically | Spenley et al ( 1993 (]. 
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3.1.2 Fracture at constant applied stress 



Another approach to study the failure of gels is to apply a constant stress and to follow 



the resulting shear rate as a function of time | Skrzeszewska et al (2010 1]. Above the 
critical fracture stress a ~ Go, fracture occurs immediately as revealed by a dramatic 
increase of the measured share rate. Below the critical fracture shear stress, a delayed 
time is observed before fracturing, a phenomenon well documented in solid materials 
IZhurkov (1965); SantuccTet al| (|2007j)] and observed also in colloidal suspensions 



[Sprakel etal(2011 



Divoux et al| (|2010|>], wormlike micelles JOlsson et aTj ( |2010) ] 



and soft elastomers | Bonn et al ( 1998) >] . The delay fracture time f/, varies significantly 



from one experiment to an other one indicating the stochastic nature of the fracture 
mechanism. With decreasing stress, the average delay time increases <%> rapidly. 
Several theoretical interpretations have been proposed, that will be reviewed in Sec- 
tion|5]l. They predict different delay fracture times. However, the intrinsic stochastic 
nature of the fracture mechanism implies a very large uncertainty of < % > that does 
not allow to conclude definitively about the relevant theory. 




10' 

time [s] 



Fig. 3 From Skrzeszewska et al 



(2010 



]: Delayed fracture of a transient polymer network formed by 
telechelic polypeptides, for different applied stresses. After a certain waiting time which varies from one 
experiment to another, the gel breaks, leading to a very rapid increase of the overall strain. 



3.2 Fractures in Hele-Shaw cells 



Fracture-like flow instabilities that arise when a fluid is injected into a Hele-Shaw 
cell filled with an aqueous solution of telechelic polymers has been reported by sev- 
eral authors (Zhao and Maher|dl993); llgnes-Mullol et all dl995b; IVlad et all (11999); 
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(b) 





Fig. 4 From Zhao and Maher 1993 ]: Four typical patterns of a 2.5% solution of hexadecyl end-capped 
polymers of molecular weight 50 700 at injection rates (a)Q = 1 .0 mL/min, (b) Q = 1 .OmL/min, (c) Q = 5.0 
mL/min, and (d) Q = 20 mL/min. The crossover from the viscous-fingering pattern to the fracturing pattern 
is seen in (a) and (b). 



Mora and Manna (2010)]. Zhao and Maher |Zhao and Maher ( 1993)] showed using 
a radial Hele-Shaw cell, that for transient networks of associative polymers, there ex- 
ists a threshold injection rate, below which the pattern instability is typically viscous 



fingering (Saffman-Taylor instability) [Bensimon et al (1986)], and beyond which, 
the patterns resemble fracture patterns observed in brittle materials (Figure |4|. Note 
that analogous fracture patterns have been also observed in other classes of com- 
plex fluids like clay suspensions | Lemaire et al| ( [T99T] )] or lyotropic lamellar phases 
IGreffier et al ( 1998 1]. The difference between a fingering pattern and a fracture pat- 



tern is drastic and can be quantified by calculating the mass fractal dimension of a 
given pattern: it drops drastically from ~ 1.70 to ~ 1.0 near the transition threshold. 
Interestingly, the fingering/fracture transition is not observed for the corresponding 
entangled homopolymers solutions, indicating that the fingering-fracturing transition 
is a direct manifestation of an associating-network effect. A characteristic, Deborah 
number D e = Xt hin /Xf — (xQ) / '(b 2 L ceU ) can be defined, where — 

l/Ytiun is the 

inverse of the shear thinning rate of the transient network and not its relaxation time 
T, Q is the injection rate, b the thickness of the cell and L ce u, some characteristic 
length scale of the radial cell. The transition from fingering to fracture-like behavior 
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Fig. 5 From Ignes-Mullol et al 1995 ]: Fracture involving an associative polymer polymer of molecular 
weight 17400. Top: viscous finger at early times. The shape is that of a Saffman- Taylor finger. Middle: 
the fracture-like pattern emerges cleanly from the unperturbed finger. Bottom: tip velocity as a function of 
tip position during the growth). 



appeared roughly at the same Deborah number under variation of polymer molecular 
weight or polymer concentration. 

Dynamics fracture-like flows instabilities of associating polymers solutions has 
been also characterized in a rectangular Hell-shaw cell by measuring the crack tip 
velocity JIgnes-Mullol et at ( 1995j>] or the tip mobility which relates the velocity of 
the tip to the pressure gradient | Vlad et al (1 999)]. The resulting pattern of the flow 
instability is highly dependent on the concentration, molecular weight or architecture 
of the associative polymer solutions. For low molecular weight polymer, the onset of 
fracture-like pattern evolution is accompanied by an abrupt jump in the tip velocity 
followed by a constant low acceleration (Figure [3]). At high molecular weight, the 
transition is blurred by a meandering regime in which the crack tip meanders from 
side to side and fluctuates with a characteristic frequency proportional proportional 
to the fluid injection rate. 

Mora and Manna [Mora and Manna| ( |2009l|2010||2012| l] revisited both theoreti- 
cally and experimentally the Saffman-Taylor instability for a model viscoelastic fluid, 
i.e. a upper-convected Maxwell fluid, a versatile experimental realization of which 
consisting a transient network of oil in water microemulsion droplets droplets re- 
versibly linke d by telechelic polymers [Fil ali et al| (fl999)]. They prove [M ora and| 

b 2 (-VPfl 2 



Manna 



ness of 



(2010 1] that a unique dimensionless parameter f = G ri j 2 (b is the thick- 
mere Hele-Shaw cell, VP is the applied pressure gradient, F, the surface 
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Fig. 6 From | Mora 201 If ] : Aspiration of a Maxwell fluid in a radial Hele-Shaw cell: (a) Elastic pattern 
formation for i = 5, (b) Fracture-like pattern for X > 10. 



tension of the network and Go is the zero shear modulus), controls all elastic effects 
and particularly, the fracture transition. Note that f does not depend on the fluid vis- 
cosity but on the elastic modulus. The theory describes the transition from viscous 
to elastic fingering instability and shows that there is a divergence of the maximum 
growth ratio of the instability, which occurs for a finite value of the control parame- 
ter f = 10.2. Experimentally (Figure [6]), the main results is the fracture-like pattern 
observed above a critical value xf p ~ 10 for every fluid, very close to the predicted 
theoretical value. So, the fracture is a direct consequence of the predicted blow up 
and the existence of the fracture-like pattens is a consequence of the elasticity of the 
complex fluid and takes place within the linear regime. 



3.3 Fractures in extensional geometry 



3.3.1 Filament stretching rheometry 

Extensional properties including fractures of associative polymer networks remain 
virtually unexplored until the development of filament stretching rheometry )Anna| 
etaIl([200T]);|McKinley and Sridhar| ( [2002] l; |Anna and McKinley] ( p008] l]. Tripathi et 
al. | Tripathi et al| ( |2006| l] investigated the non-linear extensional rheology of end- 
capped urethane coupled poly-(oxyethylene) (HEUR) using a filament stretching 
rheometer. A nearly cylindrical sample of transient network fills the gap between 
two rigid circular plates that are moved apart to a final separation with an exponen- 
tially increasing separation profile L(t) — Loexp(ef), where e is the strain rate, the 
evolution of the tensile force and mid-filament diameter being measured simultane- 
ously. The authors varied both the polymer concentration and the molecular weight 
of the polymer. They observed that, for sufficiently low polymer concentration (suffi- 
ciently low network strength), and at high strain rates and Deborah number (defined 
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( 2006 1] : Images showing rupture of 2.4wt % HEUR22-2, for an applied strain 



as De = re, where z is the Maxwellian relaxation time), a defect appears in the fila- 
ment and a tear propagates rapidly across the filament leading to a complete rupture 
in two distinct domains followed by an oscillatory damping elastic recoil of the two 
pieces toward the two end planes (see Figure [7). The failure of the filament occurs 
before any significant necking has occurred. This type of rupture instability has also 
been observed in polymer melts flJoshi and Denn|( 2004i] and entangled solutions 
of wormlike micelles [Rothstein (2003); Bhardwaj et al ( 2007) )] . At lower Deborah 
numbers, the fluid filament pinches off through a visco-capillary or elasto-capillary 
thinning, and occurs after a considerable necking of the filament. These two dif- 
ferent rupture scenarios (elastic rupture versus necking and pinch-off) observed in 
transient networks and more generally in polymeric liquids are reminiscent of brittle 
versus ductile fracture behaviors |Tabuteau et al (2008)] observed in solid materials, 
reviewed in section [4] 

Tripathi et al | |Tripafhi eTa l (2006)] interpret the failure behavior of the transient 
network filaments in the context of the Considere construction | Considere] ( j 1 885[ )] , 
originally developed in solid mechanics to quantitatively predict the critical Hencky 
strain to failure, despite the doubts about its applicability to the prediction of the 
inception of necking in extensional flow [Petrie (2009)]. The original Considere con- 
struction is that, if the graph of force against stretch ratio has a maximum, there is the 
possibility of two stretch ratios corresponding to the same force, which can be associ- 
ated with the formation of a neck in the solid sample, leading to the failure. Renardy 
[Renardy ( 2004[ l] noted that fluid threads described by constitutive models in which 
the extensional viscosity passes through a maximum may undergo a purely elastic 
mode of necking failure in which surface tension plays no role. The observations of 
Tripathi et al of filament necking and rupture processes observed experimentally ap- 
pear to be connected to an instability resulting from saturation in the tensile stresses 
at a critical deformation rate. It has been associated with a maximum in the exten- 
sional viscosity as a function of rate of strain in the constitutive rheological model 
used to described transient networks, so that this maximum of extensional viscosity 
is claimed to be an analogue for liquids of the Considere criterion for solids. How- 
ever, such an analogy to interpret the failure of transient networks requires great care 
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and at worst leads to confusion and misleading conclusions as developed in the re- 
view of Petrie |Petrie (2009 )], due to the use of too crude constitutive models for the 
complex fluids and a strong dependence on the details of the extensional flow. Re- 



cently, Fielding | Fielding (201 1 1] finds, contrary to the Considere criterion, the onset 
of instability to relate closely to the onset of downward curvature in the time (and so 
strain) evolution of the z component of the molecular strain, for extension along the 
z axis. 

3.3.2 Fracture of pendant drops 

The fall and failure of a viscous fluid droplet from a faucet has been extensively 
studied both experimentally an d theoretically (see the reviews by Eggers and Viller- 
maux [Eggers and Villermaux (2008)] and Eggers [Eggers (1997)]). However, for 



non-Newtonian fluids, such a surface free flow lead to highly non-linear fluid motion, 
in particular strong stretching that are usually non-well tested by Non-Newtonian 
constitutive laws, and few quantitative experiments have been available. Smolka and 
Belmonte | Smolka and Belmonte (2003 )] reported the filament dynamics and rupture 
of pendant drops of viscoelastic wormlike micellar fluids: the rupture is consecutive 
of a necking and pinching off instability followed by a partial retract of the free fila- 
ment ends, is also reminiscent of a ductile fracture behavior observed at low Deborah 
numbers in filament stretching rheometry experiments performed on transient net- 
works | |Tripathietai| ( p006l )]. 

The single report (to the best of our knowledge) of fail and rupture of pendant 
drops of transient network has been done by Tabuteau et al [Tabuteau et al (2009a)]. 
The experimental system they used is a model network consisting of oil in water mi- 
croemulsion droplets linked by telechelic polymers (described in section 2), which 
has the peculiarity to exhibit a pure linear Maxwell behavior as long as it flows. A 
drop emerges at a constant rate of 2mL/h from a plastic tube of 2.6 mm diameter. 
The drop of fluid begins to fall when its weight exceeds the surface tension retaining 
force; a filament is formed and stretched by the failing drop until a crack nucleates 
somewhere around the filament (Figure [8]) and then the fracture propagates perpen- 
dicularly to the filament. The total duration of the fall before the opening is typically 
on the order of 50 to 200 s depending on the sample, much larger than the relax- 
ation time of the fluid (~ Is), itself much larger than the total duration of the rupture 
event (~ few ms). After that the complete rupture was achieved, an elastic retraction 
of the remaining pendant piece of the filament is observed. Upon growing, the frac- 
ture exhibits a parabolic shape expected for an elastic solid breaking under tension 
HGreenwood and Johnson| ( |1981| );[Barenblatt]( |1962| )]. The authors measured the di- 
ameter of the drop where the fracture occurs and weighted the mass of the falling 
part, so they can measure the rupture stress a c as a function of the shear modulus of 
the gel, that can be varied by varying the polymer concentration. Moreover, the au- 
thors performed several measurements of the rupture stress of two samples with the 
same composition (oil mass fraction: 10% and mean number of polymer sticker per 
droplet: 6) and the same shear modulus Go = 1200 Pa but with two different lengths 
of the hydrophobic alcane chains constituting the stickers of the telechelic polymers 
(21 methyl groups and 18 methyl groups respectively). The relaxation time of the 
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Fig. 8 From Tabuteau et al 201 1} ]: Sequence of images of the fall of a drop of a transient network. 
Regime 1 (drop formation) corresponds to images (a) and (b) and regime 2 (viscous fall) corresponds to 
images (c) to (g). Image (g) corresponds to the beginning of the brittle fracture regime. At the beginning of 
regime 2, the radius of the filament is defined as R* = D* /2 and Lq* is the length of the drop. The white 
scale bar corresponds to 2.5 mm. 



C21 sample is five times larger than the relaxation time of the Cis sample. The mean 
value of the rupture stress is found to be identical for the two samples, with also the 
same relative standard deviation. This indicates that the rupture stress is independent 
of the binding energy. 

An other important point reported in [Tabuteau et al (2009a)] is that the rupture 
occurs in the linear regime a c ss Go, so that the critical rupture stress is close to shear 
modulus of the network (see Figure 10 1. The thermally activated fracture model de- 
veloped in section 5.2 allows a quantitative interpretation of these behaviors. One of 
the key points to understand this mechanism of fracture is the bond reversibility and 
the corresponding relevant ultra low interfacial energy needed to nucleate the crack. 
This interfacial tension results from the loss of conformational entropy of polymeric 



bonds near a crack interface and is typically on the order of few jiNm |Sprakel et al 



(2007)]. The origin of this ultra low interfacial tension will be discussed in section 
\5W 



3.4 Fracture propagation 



The propagation of fractures in transient networks is less documented |Ignes-Mullol 
|et al| ( |1995| l; |Vlad et al| ( |1999| ); |Tabuteau et al| ( |2011| i] than in permanent gels, such 
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Fig. 9 




T 



Pictures of the propagation of a fracture across a pendant filament 



of a transient network from the right to the left. The time left to achieve complete fracture of the filament 
corresponding to each picture labelled with a letter is: a (8.50 ms), b (5.16 ms), c (2.67 ms), d (2.00 ms), 
e (1.00 ms), f (0.33 ms), g (0.17 ms) and h (O.ms). The last picture on the right shows almost all of the 
elongated drop, with the crack being well developed. The white scale bar corresponds to 0.1 mm. 
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Fig. 10 From Tabuteau et al 2009a ] : Plot of the rupture stress coming from the pendant drop experi- 
ments (o"(,, circles) and the expected values coming from Eq.^](o~i, line) as a function of the connectivity 
r. The vertical dashed line corresponds to the percolation threshold (r=3). The inset gives the variation of 
the experimental rupture stress with the shear modulus, corresponding to the different connectivity. Note 
that the observed deviation of o~[ from its linear dependence with the shear modulus Go observed at small 
r, originates from the vicinity of the percolation threshold |Filali et al 1 1999)]. 
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as covalent hydrogels [Tanaka et al (|2000|l; Livne et al (2007 2005 1], biopolymer 
hydrogels ( gelatin [[ Baumberger et al (|2006b|a)>; |Baumberger and Ronsin ( 2010a| l] 
or alginates [Baumberger and Ronsin (2010b| |2009| )]) and permanent self-assembled 
triblock copolymer gels | Seitz et al (2009)]. If one excepts some reports on dynamics 
fracture of associating polymers in Hele-Shaw cells |lgnes-Mullol et al ( 1995} ); Vlad 



|et al| ( |1999] )] already reviewed in section 3.2 the systematic analysis of morphology 
and propagation of a crack in a transient network has been reported very recently in a 
single publication [ Tabuteau et al ( 20 1 1 )] . The fracture propagation arising in pendant 
drop experiments has been tracked by high speed velocimetry [ Tabuteau et ai| ( |20lT| )] . 
This configuration allows for an excellent reproducibility of the fracture initiation and 
propagation, as well as a pure elongational stress condition due to the lack of contact 
with solid interfaces near the fracture region. The fracture of the fluid happens in two 
steps. The fracture initiation step, reviewed in the previous paragraph and reported in 
| Tabuteau et al (2009a|] was shown to be governed by the thermally activated nucle- 
ation of a critical crack in the polymer network. The second step, consisting in the 
rapid propagation of a brittle fracture in the fluid, was analyzed in detail and shown to 
be energetically governed by the surface tension of the solvent (oil- in-water droplet 
microemulsion). Indeed, from the beginning of the propagation up to the complete 
fracture (i.e., when the sample is separated in two parts) the fracture profile exhibits 
a parabolic shape (Figure [9]) as expected for an elastic solid breaking under tension 
| Irwin ( 1957 »] . These observations were confirmed by quantitative analysis of the 
fracture profile u(x) on the crack propagation across the sample (Figure 1 1 ). 

Interestingly, the authors show that, because of the low value of the elastic mod- 
ulus of the transient network, and as in soft solids [Seitz et al (2009 1], the crack tip 
region exhibits very large deformations that require the use of finite elastic theories 
[Hui et al (2003])] for a proper description. The non-linear character of the problem 
requires the fracture energy to be estimated by the J-integral method |Rice| (l968 )] 
and the strain energy release can be related to the curvature radius of the tip R t t p 
through £f « J = nGoRtip/2. The brittle character of the fracture propagation is con- 
firmed by the fact that the measured fracture energy is independent of crack velocity 
V and equal to the Dupre work, that is the reversible work needed to form two new 
air/gel surfaces 2T S (dry fracture), without any bulk viscous dissipation. According to 
the viscoelastic trumpet model of de Gennes |de Gennes| ( |1996| )] applied to a Maxwell 
fluid [Saulnier et al (2004)], the absence of bulk viscous dissipation was justified by 
the short length L of the crack (less than the filament diameter) in relation to the 
characteristic length Vt above which the viscous dissipation starts to become effec- 
tive (Figure 12). An extension of the trumpet model to the hyper elastic case allows 
to confirm the weakness of the dissipated energy in terms of a new relevant length 
scale Rfip ~ ^ /Go, of the same order of magnitude as the observed crack lengths. 
In agreement with this interpretation, the measured crack opening profiles present a 
constant parabolic shape. 

Despite the brittle character of the propagation, the authors report a slow velocity 
of crack propagation V sa mm/s <C cr, cr « yjGo/p on the order of lm/s, being the 
speed of sound in the medium, expected for a brittle solid. The authors propose a 
microscopic model for thee scaling of the velocity V of the fracture to be governed 
by the time scales of the local crack tip debonding processes (Figure[T3]l, and namely 
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Fig. 11 From Tabuteau et al 



2011 



]: Fracture profiles u(x) for different times before the break-up in 
the fracture moving frame, corresponding to the pictures of figure [9] The black line is a parabolic fit 
corresponding to u{x) = a\f{x) with J = nGoa 1 /4 and J = 2r s . We report only the profile for L < 0. \Dq 



Fig. 12 From Tabuteau et al 




*■ X 



A schematic representation of the space and time scales associated 



to a crack of length L moving with velocity V in a Maxwell fluid of relaxation time T according to the 
trumpet model| de Gennes ( 1996 1]. A small microscopic nonlinear zone of length / is represented in black. 
The behavior of the material is solid like at scales smaller than Vr, then fluid like at larger scales, and 
accordingly the fracture profile scales differently with the distance to the crack's tip in the two different 
regions. 
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Fig. 13 From Ref Tabuteau et al 



201 1 1]: Cartoon of the viscous relaxation mechanism of a bead at the 



tip of the fracture. The polymer bridge between beads (1) and (2) just debonded, forming a loop on bead 
(2). Bead (1) thus experiences the spring-back force /(f) due to the gel under tension; this will lead to an 
increasing extra tension on bead (3) and crack propagation at velocity V. 



by the elastic relaxation time of the oil droplets after a debonding event under the 
action of the viscous drag of the solvent (Reynolds number R e < 1). This leads to an 

upper limit for the estimation of the crack velocity V < 0.25 2%^ w b 1 wnere Go is the 
shear modulus of the gel, b the radius of the droplet and % the mesh size and rjw is 
the viscosity of water. The control of the crack velocity by network/solvent friction 
has been already observed by Baumberger et al. [Baum berger et al (2006b a)] in the 
viscoplastic fracture dynamics of thermoreversible gels. 



4 Brittle vs Ductile fractures 

Ductility versus brittleness characterizes the nature of fracture mechanism in solids 
under stresses when irreversible damages take place. Brittle solids suddenly break in 
the elastic regime above a critical stress, whereas ductile solids plastically deform, 
before fracturing. Ductility implies that some specific mechanism exists, allowing 
local rearrangements and creeping while the micro-structure is on the average pre- 
served. Speaking of brittle fractures for viscoelastic fluids seems a priori paradoxical 
since, like all liquids, they accommodate arbitrarily large deformations after a finite 
relaxation time: in that sense they are infinitely ductile. However, since the pioneering 
works of Vinogradov and coworkers in the 1970s and 1980s, it has been recognized 
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that polymers melt can fail either in a ductile or a brittle mode, as discussed in details 
in | Wang and Wang (2010 1; Boukany et al (2009)]. The switch of the failure mode 
from yielding to elastic rupture is controlled by the Rouse Weissenberg number. 

The observation of fluid fractures in extensional geometry clearly reveals two 
distinct rupture scenari (elastic rupture versus necking and pinch-off) observed in 
viscoelastic fluids and reviewed in section |3.3| that are reminiscent of brittle ver- 
sus ductile fracture behaviors observed in solid material. Interestingly, two transient 
networks (one consisting of cross-linked polymers the second of entangled living 
polymers [Dreiss (2007 )], the other having the same simple Maxwell linear Theologi- 
cal behavior exhibit the two distinct different fracture behaviors as shown in pendant 
drop experiments [Tabuteau et al (2008)] (Figure 14 1. Intuitively, the fracture should 
occur when the material does not have enough time to relax, i.e. when the flow time 
scale is faster than the relaxation time. This condition is equivalent to the condition 
that the Deborah number is larger than 1 | Gladden and Belmonte (2007)]. However, 
as shown in Figure 14 and in [Tabu teau et al| ( |2008[ )], this appealing simple picture is 
a bit too naive. Bhardwaj et al |Bhardwaj et al (2007)] observed both types of frac- 
ture in solutions of entangled micelles and suggested that the smooth pinch-off mode 
observed for the more concentrated sample might be the results of structural modi- 
fications of the solution of entangled wormlike micelles. On the other hand, Sprakel 
et al | Sprakel et al ( 2009b])] investigated flow failure modes for transient polymer 
networks using particle-based simulations and showed transition from shear band- 



ing to melt fracture upon changing the polymer concentration. Tixier et al |Tixier 



|et al| ( |2010| l] reported the existence of a new class of self-assembled transient net- 
works made of surfactant micelles of tunable morphology (from spheres, to rodlike 
to wormlike, Figure 15 1 reversibly linked by telechelic polymers, designed to achieve 
a comprehensive study of the brittle to ductile transition in self-assembled networks. 
Flow curves of the transient networks of tunable morphology suggest that one can 
associate three domains as the micelles growth with distinct failures modes: a brittle 
mode, an intermediate mode and finally a ductile/shear banding mode, as the micelles 
grow. 

Coupling rheology and time-resolved structural measurements using synchrotron 
radiation [Ramos et al (2011 1], Ramos and Ligoure show that the bridged micelles 
can align under shear and proposed the emergence of strong fluctuations of the de- 
gree of alignment of the micelles as a structural probe of a fracture process under 
shear (FigurqTo). The transition from a regime where the fluctuations of the are weak 
to a regime where the fluctuations are strong is rather sharp and enables to define a 
critical shear rate for the onset of strong fluctuations Yfracture- By analogy with solid 
materials, and mapping the strain in solids to the shear rate of viscoelastic fluids under 
shear, the authors therefore propose to compare the shear rate for the onset of fracture 
Yfracture and the shear rate for the onset of non-linearity Ythin ■ They define A as the 
normalized difference between the two critical shear rates A = (Yfracture — Ythin) /Ythin- 
This parameter A is analogous o the amount of plastic deformation accumulated 
before fracture in ductile materials. So, the critical shear rate for fracture is either 
very close to the shear rate at which the sample departs from its linear behavior (for 
"brittle-like" samples) or significantly larger than the shear rate at which the sam- 
ple departs from its linear behavior (for "ductile-like") samples, and the ability of 
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Fig. 14 From Ref Tixier et al 2010 ]: Middle: Shear stress relaxation function G(f) for two model ex- 
perimental systems of transient networks (bridged microemulsions and wormlike micelles). Both systems 
exhibit the same Maxwell behavior. Left: Fall of a droplet of the bridged microemulsion sample under the 
influence of gravity. Photographs taken at times (0, 10 ms, 20 ms) from the left to the right. The diameter 
of the filament is ss 5 mm. Right: Fall of a droplet of a wormlike micelles sample under the influence 
of gravity. Photographs taken at times (0, 66 ms, 67 ms) from the left to the right. The diameter of the 
filament is Rs2 mm. 
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Fig. 15 From Ref Tixier et al 



2010 



: Schematic view of the structure of transient networks with nodes 
of tunable morphology. Upon variation of the co-surfactant to surfactant molar ratio of the micelles, the 
morphology of the nodes evolves from spherical (left) to rodlike (middle) to wormlike (right). 



transient networks to exhibit a ductile-like behavior should be related to its ability to 
rearrange its structure at large scale under flow, before relaxing its elastic energy. 



5 Theoretical approaches of crack nucleation in transient networks 



Experiments reported by different groups on various transient networks but with sim- 
ilar structures and microscopic dynamical properties lead to several interpretations 
that do not seem at first glance compatible, leading to some confusion | Skrzeszewska| 
etallpOTOfrrBerret and Serero| ( |20^ ; |Tabuteau et al| p008] |2009a[ |20 1 1) ]. Sections 



5.2 |5.3 and |5.4| are devoted to theories for the nucleation of brittle fractures, and to 



the confrontation of these theories with experimental results. Section 5.2 deals with 



the thermally activated crack nucleation model (TAC model) [Golubovic and Feng 
( |1991| l; |Pomeau ( 1992} ]. In section 5.3 the activated bond rupture model (ABR 



model) pvans and Ritchie| ( [T997] l; |Chaudhuryl ( [T999| ; |Skrzeszewska eTal] ( [20T0l i] is 
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: (a), (b): Shear stress (open symbols) and Contrast (full symbols) 
of the applied shear rate, y, in the steady state. The two sets of data in (a) and (b) correspond to two 
independent measurements. The bars for the contrast data correspond to the standard deviation. (c),(d): 
Time evolution of the contrast in the steady state as a function of y. (e),(f) Normalized fluctuations of 
the contrast (black symbols) and viscosity (gray symbols) as a function of y. (a), (c), (e): brittle fracture 
behavior; (b), (d), (f) ductile fracture behavior. 



introduced. The starting point of these two models is the Griffith theory for fracture 
in brittle materials [Griffith ( 1921 1], where a critical length for crack growth is ex- 
pressed as a function of the applied load (Section |5T| i. This length is reached by the 
thermal fluctuations in the TAC model, whereas it results of an acceleration in the rup- 
ture rate of the bonds near a crack tip in the ABR model. A third model, that we call 
Self healing and activated bond rupture nucleation model (SH-ABR model) ]Mora| 
( |201 l| l] is introduced in section [54] It is based on the balance between cross-links 
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disconnection and reconnection [Mora ( |2011| l] to control crack nucleation. Finally, 



the validity of these models is compared and discussed in section 5.5 



5.1 Griffith Fracture Criterion 



The Griffith's theory for fracture in brittle materials | Griffith ( 1921| l] starts up with 
preexisting micro-cracks and considers the energy required for them to grow sponta- 
neously under a given constant stress. An interfacial energy has to be paid to break 
more cohesive bonds, counterbalanced by the bulk elastic energy released by the 
opening of the crack. For a given stress a the Griffith free energy cost W for a crack 
of size L reads, in 3D: 



W(a,L) = -L 2 F s -a^- : 



(1) 



with L the size of a disc-like crack, F s the cohesive energy per unit area, Go is the 
shear modulus and a ~ 1 a constant depending on geometrical factors. A typical 



example of the Griffith potential as a function of L is shown on Figure 17 The stress 
a being fixed, W(o,L) reaches a maximum: 



W max {a) 



6tz F?Gl 



a- 



(2) 



for a crack size (the so called Griffith length) 




(diameter of the crack) 



Fig. 17 Sketch of the Griffith's free energy of a crack as a function of the crack length L. Note that this 
curve can be scanned in both directions (increasing or decreasing L) for self-healing materials, in contrast 
with usual non-self-healing materials. 



L c (a) = w-. 

a o 2 



(3) 
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For a crack size larger than L c ,dW jdL < 0, leading to a catastrophic growth of the 
crack (fracture), in contrast with a crack size smaller than L c . Thus, a stressed sam- 
ple appears to be in a metastable state as long as no crack with a size larger than L c 
nucleates [Golubovic and Feng ( 1991 1; Pomeau (2002)]. 



5.2 Thermally activated crack nucleation model (TAC model) 



A material submitted to a constant stress breaks after a given waiting time, that de- 
pends on the applied stress. Delayed fractures have been reported for a large variety 
of materials, ranging from 2d-microcrystals |Pauchard and Meunier (1993 1998|l] 



to heterogeneous materials [Guarino et al (1999i; Ciliberto et al (2001); Guarino 
et al ( |2002| i] and gels [Bonn et al| ( | 1998 1]. More recently, experimental evidences of 



delayed fractures have been reported in transient networks |Tabuteau et al (2009a); 
Skrzeszewska et al] (|2010|l] and a good agreement has been found with predictions 



of the TAC model [ |Pomeau| ( [1992 )] , whose starting point is the Griffith theory. Sev- 
eral auth ors ||Golubovic and Feng| ( |1991| l; |Pomeau| ( [l^92] i; |Buchel and Sethna] ( |1997| l; 
|Pomeau| ( |2002[ l] have proposed that the Griffith's energy (Eq.[T]i is an energy barrier 
of height W max ((j). Then, the waiting fi should follow the Kramer's theory [Hanggi 
[etaIl ( [T990l l]: 



(h) °= exp 



k B T 



(4) 



where (.) is for averaging over the stochastic distribution of fracture waiting times. 
The choice of an energy barrier implicitly assumes that a crack can reversibly explore 
states with different crack lengths between the initial and the critical ones. 

A key characteristic of crack nucleation in transient networks is the intrinsic re- 
versibility of the process. This directly comes from the transient nature of the junc- 
tions: new junctions can randomly appear between the two opposite sides of a crack. 
This is in fundamental contrast with what happens in non-self-healing materials, 
where the crack length irreversibly increases, making Eq.|4]irrelevant for solid materi- 
als. Recent models can overcome this specific difficulty of non-self-healing materials 
ISantucci et al1 ( |2003][2007l >]. 

A direct way to check predictions of the TA model is to measure the average 
fracture waiting time for different values the applied stress. Putting Eq. [2] in Eq. |4] 
one obtains: 

exp (^) 4 , (5) 

where the prefactor Ti varies much more gently as a function of a, than the expo- 
nential for a G\. The variation of the fracture waiting time given by Eq.|5]is very 
stiff when the stress is close to the characteristic stress G\\ the waiting time becomes 
very long as soon as the stress decreases below <Ti, and becomes extremely short if 
the stress increases beyond <7i . This sharp variation makes difficult the experimental 
determination of the waiting time variations with respect to the applied stress: the av- 
erage waiting time should vary over several decades to check Eq. [5] Skrzeszewska et 
al. | |Skrzeszewska et al| ( |2010| l] have measured the fracture waiting time for transient 
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polymer networks formed by telechelic polypeptides, with nodes made by the self as- 



sociation of three collagen-like triple helices [Skrzeszewska et al (2009 1]. The sample 
is put in a rheometer equipped with a Couette geometry, and a given constant stress 
is applied. The measured waiting time is defined as the elapsed duration between 
application of the stress and the instant when the measured shear rate jumps sharply 
by several decades (Figure|3j. An important feature emerging from these experiments 
is the great dispersion for the measured values of the waiting times. This dispersion 
is inherent to the stochastic feature of the subcritical fracture, tf, is plotted in Fig- 
ure 18 for different values the imposed constant stress. According to available data, 



the standard deviation reached one decade. % is plotted as a function of (a /Go) (in 
semi-logscales, Figure [T8)-Left, Go being the shear modulus). Note that the data point 
corresponding to /Go = 3 comes from a strain versus time measurement whose ap- 
pearance is significantly different from the others (smallest applied stress of Figure[3]l, 
which is a possible manifestation of the phenomenon of wall slipping at high stresses 
| de Gennes] ( |2007| l; Erk et al| ( |2012[ )]. However, even excluding this data point, the 
agreement between the experimental data and the prediction is not sufficiently clear 
to state that the TAC model provides an explanation for the existence of the observed 
delayed fracture as explained in section [33] 



i * 





0.2 0.3 

WGo) 




a/G 



Fig. 18 Data from Skrzeszewska et al (2010 ]: Measured fracture waiting time (tt) as a function of the ap- 
plied stress. Left: t/, as a function of (Go/ct) 4 . The straight line represents the best fit for <oe"P /°") 4 ) 

(o"[ /Go = 1.74 ± 0.08). Right: representation in a semi-logarithmic scale. The straight line shows the best 
fit for an exponential decay to exp(— cr/a?) (Gi/Gq = 3.3 ±0.6). 



Another way to check a basic prediction coming from the TAC-model is to mea- 
sure the characteristic value <7i of the stress and to compare it to the expression fol- 
lowing the model. From Eqs.(|2]|5]l <J\ reads: 

^(^)" >c » 2 



<7i can be interpreted as follows: the fracture propagation is subcritical for a load 
smaller than <3\. Otherwise, the fracture nucleates instantaneously. Oi can be esti- 
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mated as follows. Let S, be the typical mesh size of a transient network. Follow- 
ing Green and Tobolsky [Green and Tobolsky (1946)], the shear modulus is Go — 
akgT /^ 3 , where a, of order unity, depends on the network topology. In order to esti- 
mate the value of F s , we consider a micro-crack of the network. First of all it is im- 
portant to note that micro-cracks in these networks are completely wet by the solvent 
and thus only the polymer network contributes to F s (see Figure 19 1. The energy cost 
for creating a network-solvent interface originates from the depletion of the chains 
monomers at the crack interface. We use the well known analogy between transient 
networks and polymers solution in th e semi-dilute regime : the cost in free energy per 

where a', of 



unit surface due to the depletion is [Joanny et al 
order unity, also depends on the network topology. 



1979i] F s 



U £2 i 



hen, Eq. |6]yields: 



/ 6na 



'V 



V a 2 



1/4 



k B T 

£3 



/ 6na'- 



1/4 



G . 



(7) 



The value of the prefactor being not far from unity, it appears that <7i is close to 
the shear modulus. Particle-based simulations where 0\ is plotted versus polymer 
concentration corroborate Eq. [7] ISprakel et al (2009b)]. Note that in Eq. |7J ffi is 
independent of the relaxation time and of the viscosity. It is therefore independent of 
the shear rate. 



The best fit (Figure 18 1 of the data of Skrezesezewska et al. | Skrzeszewska et al 



( 2010| l] with Eq. 13] (excluding, as explained above the data point corresponding to 
O = 3Go), gives ^ ~ 1 .74 ± 0.08, in agreement with Eq. |7j 

Tabuteau et al. have measured <3\ using pendant drop experiments, as introduced 
in section [3. 3. 2| The authors used a model system made of oil-in-water microemul- 



sion droplets connected to each other by telechelic polymers (see Figure 19 »| Filali 



et al ( 1999)] at different concentrations in order to vary both the relaxation time and 



the elastic modulus. In pendant drop experiments, the increase in stress with time is 
large enough to neglect the probability of breaking due to a subcritical fracture stress. 
Therefore the characteristic stress <7i in Eq.|7]is equal to the stress at the detachment 
in pendant drop experiments. The structure of this transient network being well char- 
acterized (particularly the network connectivity), it is possible to obtained an explicit 
expression for coefficients a and a' in Eq. [7] | Tabuteau et al ( 2009a i] for this spe- 
cific material. The agreement between measured critical stresses and the theoretical 
prediction is satisfactory (Figure 10 1. 

This result confirms the physical basis of the interfacial cohesive energy and al- 
lows the authors to conclude that the observed fractures are brittle. Note that the 
calculated value for the Griffith's length is, for all the samples, of the order of three 
times the mean distance between two oil droplets. This is consistent with the Griffith's 
theory. Furthermore, the measured critical stress has been found to be independent of 
the length of the stickers and therefore of the junction energy (note that the junction 
energy and the interfacial cohesive energy are here independent from each other). It is 



therefore independent of the transient network relaxation time | Tanaka and Edwards 



(1992e)]. This is consistent with the prediction of the the TAC-model for Cf\, where 
neither the junction energy nor the relaxation time do appear (see Eq.[7]i. 

This last point is in contrast with the experiments of Skrezesezewska et al [Skrzeszewska 
et al (2010)] (see for instance Figure [20]), showing that the TAC model fails -at least 
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Fig. 19 Schematic of a bridged microemulsion. The telechelic polymers can either link two oil droplets or 
loop on a single one. (Left) Before the crack nucleation (bold dashed line) polymers can bridge oil droplets 
on both side of the bold dashed line. (Right) When the crack occurs, the same polymers cannot cross the 
bold dashed line anymore and form bridges in the other directions. 



in part- for the experimental system used by these authors. This led them to focus on 
an alternative model, based on a more microscopic description of the phenomenon of 
fracture nucleation: the activated bond rupture model (ABR model). 



5.3 Activated bond rupture model (ABR model) 



Another way to describe fracture nucleation in transient networks is based on acti- 
vated bond rupture [Evans and Ritchie ( 1997| l; Chaudhury ( 1999 1; Skrzeszewska et al 
( |20T0l >]. 

In the systems we have in mind, the polymers are not entangled, so the single 
relevant time scale comes from the debonding kinetics. The free energy variation A&o 
characterizing a network junction is assumed to depend on the solvophobic length 
and the solvent quality. There is, by thermal fluctuations, a finite probability that a 
junction acquires sufficient energy to overcome the activation barrier A&o and detach 
spontaneously. Following the model of Green and Tobolsky JGreen and Tobolsky" 
( 1946)], the exit rate (l/lb) is estimated as the product of a natural thermal vibration 
frequency, V (~ 10 10 — 10 12 Hz) with exp (—A^o/ksT). Hence, one expects 



— = vexp - 
T V k sT 



(8) 



From the Eyring kinetic theory of fracture |Eyring ( 11936) ; Kausch ( 1978) )], Tanaka 
and Edwards [Tanaka and Edwards ( 1992b)] postulated that the lifetime of a junction 
also depends on the force carried by a link, whose stored energy can supply, through 
mechanical work, a portion of the activation energy. The junction lifetime T becomes: 



1 ({A%-f8) 



To exp 



fS_\ 
k B T) 



(9) 



where / is the force carried by the link under tension, and 8 a typical length on which 
the force works (for instance 8 is the length of the solvophobic ends of telechelic 
polymers). A direct implication of this assumption is that a stretched bond breaks 
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more quickly than the non-stretched ones. Tanaka and Edwards predict the mechani- 
cal behavior of such transient networks under different conditions (strain oscillatory, 
transient, constant shear rate, etc.) [Tanaka and Edwards (1992b c d)]; they also pre- 
dict non-linear effects (such as shear thinning) observed experimentally | |Tam et aT 



1998) ); Ma and Cooper ( 2001) >; Sprakel et al (2008)]. These non-linearities come 
from the reduction of bond lifetime under load (Eq.[9jl. However, because this theory 
is itself a mean-field theory with affine deformations of the network junctions, the 
existence of a fracture completely escapes this model. 



5.3.1 Constant load 

If the length of a pre-existing crack is larger than the Griffith's length, then the frac- 
ture develops instantaneously, without delay. We are considering here, an opposite 
situation, where the crack size is smaller than the Griffith's length. The crack then 
grows slowly and can be described by the model of the subcritical crack growth 



ISantucci et al (2007 1]. Because of the presence of heterogeneities due to thermal 
fluctuations in the network, the stress is inhomogeneous in the material, and then the 
force carried by a bond is not homogeneous. It is stronger where the bond density is 
low. Since the stress is concentrated near the crack tips, cross-links located just near 
a given tip will experience a higher force and then a higher dissociation rate. The ki- 
netics of breakage is accelerated, which will result in a shorter time for the size of the 
micro-crack to reach the Griffith length. Once the crack reaches the Griffith length a 
catastrophic fracture follows. As a first approximation (i.e. neglecting the creation of 
new bonds), the average time (fi) required for the micro-crack to reach the Griffith 
length has to decrease exponentially with the applied stress | Brenner ( 1962|; |Zhurkov| 
( (1955) ]: 

(r 2 )~exp^-|-^, (10) 

where the characteristic stress 02 incorporates the factor 5/(^7') but also takes 
into account the stress concentration near the preexisting crack. Note that the co- 
variation of the Griffith length (Eq.[3]) has been neglected because it varies much more 
slowly than the exponential. 



Such dependence of the fracture waiting time (Eq. 10 1 as a function of the ap 



plied stress has been observed in several types of colloidal gels |Divoux et al|(f2010 ); 



I Sprake l et al (2011)]. Skrzeszewska et al. JSkrzeszewska et al| ( |2010 i] have checked 



the validity of this approach for transient networks subjected to a constant load. They 



found a good agreement between Eq. 10 and their experimental data, as shown in Fig- 
ure 18 Right. However, it is difficult to conclude from these data what is the relevant 
model among the TAC-model and the activated bond rupture model (ABR model); 
Indeed, none of the two plots of Figures [T8]cleaiTy prevails over the other. 

Furthermore, the range of validity of Eq.[l0]is limited, since cross-links re-formation 
are not taken into account. It can therefore be so inferred that this model cannot be 
applied for small loads, typically when the fracture waiting time (t2) exceeds the 
characteristic time of cross-links re-formations. In addition, Eq. [lOlalso predicts a 
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smooth variation without any abrupt jump for the fracture waiting time regardless 



of the value of the applied stress, contrary to other models (see sections 5.2 and 5.4 



and Figure 23 1. This is in opposition with the existence of a experimentally observed 



fracture threshold (an unstressed transient network does not spontaneously break). 



5.3.2 Constant shear rate 



Sections 5.2 and 5.3.1 deal with crack nucleation for transient networks subjected to 



constant loads. In this case, the sample flows before breaking because the fracture 
does not appear instantaneously and because the material is a fluid. In this situation, 
the shear rate varies with time because it is subject to fluctuations even when the 
steady flow is still established. 

Very commonly, a shear rate is imposed in rheological tests. Some studies have fo- 
cused on fracture nucleation in a transient network subjected to a constant shear 



rate (see e.g. section 3.1.1 1 |Berret and Serero (2001 1; Skrzeszewska et al ( 2010 1]. 



Skrzeszewska et al. [Skrzeszewska et al (2010)] observed that the failure stress and 
the waiting time (before the fracture occurs) depend on the imposed shear rate (Fig- 
ure 20 1: the waiting time is found to decrease and the critical stress to increase as the 



constant shear rate is higher. Interestingly, the critical strain exhibits a non monotonic 



behavior on the share rate (Figure 20 1), very similar as what is observed in uniaxial 



extension of polymeric liquids |Malkin and Petrie ( 1997| l]. Particle-based simulations 
also show a weak but significant increase of the fracture stress with applied overall 
shear rate [Spr akel et al| ( |2009b| l]. The authors of [Skrzeszew ska et al| ( |2010| l] have 
used the ABR model to explain these results. Within a mean field approximation, the 
deformation is affine and a bridging chain in the network is stretched at a speed £, y, 
where y is the (macroscopic) applied shear rate and S, the typical distance between 
junctions. Supposing that bridging chains act as Gaussian chains with a spring con- 
stant kgT /E, 2 , the force carried by a cross-link increases as /(f) = K^fft, where 
the factor K accounts for the local stress enhancement in the vicinity of the crack tip 
| Skrzeszewska et al ( 20 1 1] . With Eq.|9]one obtains, for the survival probability of a 
bond P(t) ]Evans arid Ritchie| ( |1997| ]: 

f f , i 1 [Kk B Tjt'8 
P{t)=exp< - / df' — exp' 



The average lifetime of a cross-link is 



T ) = / P(f)df. 



(ID 



(12) 



In the high shear-rate limit, yr 3> ^§ and yr 3> 1, re-formation of cross-links can be 
simplifies and it follows that a fracture occurs after a time t' 2 ~ (tQ 



neglected. Eq. 



12 



z KSy 

and the critical strain reads Y 2 = Vi — 



In 



V 4 / 



hardening [Tabuteau et al (2009b i; Serero et al (2000 1], the stress at the rupture (o' 2 ) 



(13) 



Neglecting a possible strain 



30 



Christian LIGOURE, Serge MORA 



can be estimated as a ~ Go 7. The experimental critical strain Y 2 is plotted in Figure 



20 A for three different concentrations of a transient polymer networks formed by 



telechelic polypeptides (from [ Skrzeszewska et al ( 2010| l]). In order to make a more 
quantitative comparison between their experiments and the model, the authors have 
plotted dimensionless values of and a 2 {y^S/E? and a 2 8 / (Gq^)) as a function of 
In (yr 5 /<!;), where 8 is the estimated length of a junction and % is estimated using 
Go — kgT I '4 3 . The curves superimpose with a nice agreement with theoretical pre- 
dictions of Eq. [12] in the high shear limit. This gives a strong argument in favor of 
the relevance of the ABR model for this system, whereas the TAC-model predicts a 
shear rate independent failure stress. However, these predictions does not seem ap- 
plicable to other transient networks, where a shear stress independent failure stress 
has been reported [Tabuteau et al (2009a)]. As explained below, the main weakness 
of the ABR model is that cross-links re-formations are neglected. 
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Fig. 20 From Skrzeszewska et al 2010 ]: A- Critical strain as a function of shear rate for transient 
polymer networks formed by telechelic polypeptides at different concentrations. Lines corresponds to 
the predictions of the ABR model. (B) Renormalized critical strain as a function of the logarithm of a 
renormalized shear rate, (C) Same for critical stress. 



5.4 Self healing and activated bond rupture nucleation (SH-ABR model) 



The fracture waiting time predicted by the ABR model leads to a finite waiting time 
whatever is the applied stress, even if it is zero. Following this prediction, an un- 
stressed transient network should spontaneously break after a finite time. This in- 
consistency arises because the creation of new bonds (which could replace the bro- 



ken bonds) is not taken into account. Recently, Mora |Mora (2011 1] has proposed 



a model, the Self healing and activated bond rupture nucleation model (SH-ABR 
model), which aims to describe fracture nucleation in transient networks considering 
both activated bond rupture and cross-links re-formation to overcome the intrinsic 
limitation of the ABR model. 

It is based on the fiber bundle model, usually used to study fracture in heteroge- 
neous materials | |Coleman| ( [7958] ); |Moreno~et"all ( fT999) l; |Kun et al| pOOOl ; ICiliberto] 
|et al| ( |200T| ); |Pradhan et al ( |2010 l]: a set of fibers (or bonds) is located on a support- 
ing lattice and a random stress threshold is assigned to its elements. The set is loaded 
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and a bond breaks when its load exceeds a threshold value. In the equal load sharing 
hypothesis (or global), which is the simplest scheme one can suppose, it is assumed 
that the load carried by failed elements is equally distributed among the surviving el- 
ements of the system. The SH-ABR model deals with systems without entanglement. 
The relaxation time of a link is then assumed to be far shorter than the bonds life- 
time. In this model, a bond just once broken becomes active again: It can either take 
back the place it had or bridge another junction. Mora | |Mora| ( |201 l| l] first consid- 
ered a one dimensional model (see Figure 21 1 subjected at time / = to an external 
stress a. It is assumed that (i) following the Eyring kinetic theory of fracture |Eyring 



(|1936| l; Kausch ( |1978 1; Tanaka and Edwards ( 1992a l], the lifetime T is given by Eq. 
[9] and (ii) one end of a bond that has just escaped will reconnect over a much shorter 
duration than the (load depending) lifetime of a bond (as in most of experimental 
systems). 




Fig. 21 From Mora 2011 ]: Schematic of the one-dimensional model. The system is composed of N 
layers. The layer number ;' is connected to the layer number (' + 1 by n, bonds. The total number of bonds 
and layers is fixed once and for all. All bonds are identical. They can connect only two adjacent layers. 
One end of a bond can migrate from one layer to another. A force applied to the ends is transmitted layer 
by layer. Within a layer, this force is split evenly between the springs. 



An analytical critical magnitude of the applied stress is predicted in ID: 



03 



n Q ^k B T 



(14) 
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where no is the number of network junctions per unit volume and S, the mesh size 
of the network (the distance between two layers in the one-dimensional model, see 
Figure 21 1. Beyond this threshold, any concentration inhomogeneity in the bonds dis- 



tribution is amplified, leading to the formation and the catastrophic and irreversible 
growth of cracks. On the contrary, any inhomogeneity in the bonds distribution de- 
creases in time when the applied stress is below this threshold value. Note that the 



load dependence of the bond lifetime (formally: 8 ^ in Eq. 14 1 is here a necessary 
condition for the existence of a fracture stress threshold. 

The time (t^) at which the concentration of bonds vanishes somewhere in the one- 
dimensional system is by definition the failure time. This time is found to depend on 
the ratio " a ^ T (Figure 



22 1. It tends to infinity as 



np^ B T 



aS 



tends to one i.e. when 



approaching the failure threshold (a = 03). Of course, the failure time is not defined 

for n -^f- > 1. 
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Fig. 22 From Mora 



2011 



to infinity on approaching 



: Reduced failure time ^3 /tq as a function of yo = *°^ r . This time diverges 
le rupture threshold. Insert: Same curve in logarithmic scales. 



A similar result was found when the strain is imposed. The threshold value of the 
deformation beyond which rupture occurs is simply equal to the ratio of the critical 
stress (Eq. 14 1 with the elastic modulus of the system. The fracture waiting time is 
finite when the imposed strain is above this threshold. It decreases gradually as the 
imposed strain increases. 

Numerical simulations are used to validate the results, at least semi-quantitatively, 
for two-dimensional networks |Mora (201 1 )]. Note that the fractures under a constant 
shear rate have so far not been studied within this model. 

In conclusion, the introduction of a self-healing process in the ABR model leads 
to a true fracture stress threshold. 
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5.5 Ranges of validity and discussion 

We have presented three models that lead to different predictions for the fracture 



waiting time and for the critical failure stress (see Figure 23 1. 

When the applied stress is large, the average bond lifetime is very short (see Eq. 
|9|l and the effects of cross-links re-formation become negligible. Here, large means 



exp(— fS/(ksT)) ^ 1, i.e. o ^> 03 (where 03 is given by Eq. 14 1. Figure 23 shows 
that for a ^> 03 the decay of the fracture waiting time as a function of the stress is 
exponential, in accordance with the prediction of the ABR model, i.e. (t^) ~ (t2). So 
in this limit, the ABR and the SH-ABR models lead to the same predictions. 

In the opposite limit, where exp(— f8 / (kgT)) < 1 or equivalently a < 03, the 
effects coming from the load dependence of the cross-links lifetime are negligible 
(formally, 8 = in Eq. [§). The ABR model is no longer valid because the cross- 
links re-formation events are important. In this case the SH-ABR model predicts the 
absence of fracture (because a < 03). Within this limit, the TAC-model is the only 
one able to explain the origin of the experimentally observed fractures. 

Because (to,) — ^ °o (resp. 0) for a — > o7 (resp. °°) and (fi) — > °° (resp. Ti) for 
(7 — s- (resp. °°), (fi) is smaller (resp. larger) than (tj) for low (resp. high) stresses 



(Figure 23 1 



Let us suppose that the relevant model is the one which predicts the shortest 
fracture waiting time for a given value of the stress. If (7i < 03, then (fi) intersects 



?3) for a stress larger than Ci (Figure 23 right). As a result, (t?,) will be the relevant 



fracture waiting time for stresses for which (fi) is of order Ti (a molecular time, 
Eq. [5]), in contradiction with the implicit hypothesis of a finite breaking time. In the 
opposite case, where di > 03 (Figure [23[left), (fi) intersects (t$) for a stress smaller 
than (7i for which (f 1 ) is a priori far larger than Ti . The SH-ABR model then predicts 
a macroscopic fracture waiting time smaller than (t\). It is so the relevant model in 
this case. 

The TAC and SH-ABR models described above must be thought of as asymptotic 
models. For intermediate stresses a more complete description must take into account 
all the ingredients of both models to give a more refined transition between the TAC- 
model and the SH-ABR model. 

The condition for which the SH-ABR model is likely to adequately describe the 
fracture waiting time is (i) <7i > (73 and (ii) the applied stress is greater than the 
crossover stress for (fi) and (t^). From Eq.|7]and 14 condition (i) yields: 



feA 2 V /4 n^k B T 4 

<->0 > g ~ -jtro- (15) 



Eq. 15 leads to the approximate condition (omitting the prefactor of order unity): 

8>t;. (16) 
The TAC-model applies at any stress when (73 > Oi, i.e. for £, ^> 8. 



Let us come back to the two main experimental systems introduced above. In 
the experiments of Tabuteau et al. [Tabuteau et al (2008, 2009a)], the mean distance 
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Fig. 23 Fracture waiting time ti, as a function of the applied stress calculated for the three models. Top: 
o"i > CT3 . In this case, the curves giving (f i ) and (tj ) intersect and the corresponding waiting time is not 
molecular {tt 3> Ti). Bottom: o~3 > o~[ . (t\ ) and (13) intersect at a molecular time (Ti in Eq.|5j. 



between junctions is of the order of £ ~ 130 nm whereas the length of a hydropho- 
bic sticker is about 8 ~ 10 mm. Thus, 8/£, < 1 is consistent with the experimental 
measurements in agreement with the TAC-model. 

Because the fine structure of the triple-helical junction in transient polymer net- 
works formed by telechelic polypeptides as studied by Skrzeszewska is not well 
known, the estimation of 8/S, is less accurate. The authors took 8 ~ 9 nm, lead- 
ing to a ratio 8j% ranging from 0.18 to 0.5, depending of protein concentration. This 
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limit corresponds to the situation where it is expected that the different models should 
be coupled. However, the experimental results obtained by the authors lead in favor 
of the ABR (or more SH-ABR) model. This suggests that the ABR model is prepon- 
derant even for values 8 /E, slightly below 1 . 

In summary, we have introduced three models, each of them predicting a fracture 
waiting time as a function of the applied load. We have assumed that the relevant 
model is the model predicting the shortest waiting time. Then, it is possible to predict 
the variation of the waiting time. The threshold values for the applicability of these 
models depend on (i) the value of the applied stress, and on (ii) the ratio S/£. 



6 Conclusion and Perspectives 

Transient networks appear as simple model systems for the physicist to achieve a 
fundamental understanding of fracture phenomenons in complex fluids because of the 
unique transient character of the polymeric bonds; In this review we have focused on 
the specificities arising from the fluid and self-healing characters of these materials. 
Fractures of transient networks have been reported in several configurations: shear, 
extensional or hybrid geometries. The main results include: 

• Experimental evidences of the occurrence of true fractures, that can be unam- 
biguously distinguished from hydrodynamic instabilities in various geometries. 

• Measurements of failure thresholds, expressed as the applied stress or the ap- 
plied deformation (e.g. shear rate), and measurements of the fracture waiting time 
for subcritical stresses or deformations. Fracture waiting time as a function of the ap- 
plied stress and threshold fracture stresses have been measured by different groups. 
These measurements have led to different conclusions in apparent contradiction. Dif- 
ferent and somewhat contradictory theoretical approaches developed to understand 
the experiments have been reviewed. Two are emerging among them: the Thermally 
activated crack nucleation model (TAC model) and the Self healing and activated 
bond rupture nucleation model (SH-ABR model). Each model has a limited range of 
validity. Two characteristic lengths play a key role to discriminate between them: the 
average distance between two junctions of the network (<!;), and the size of the adher- 
ing part of the junctions (8). For £, 3> 8, the TAC model is always relevant, whereas 
for E, 8, the TAC-model is relevant for low applied stresses and the SH-ABR model 
is relevant for high stresses. This should explain the diversity of the experimental be- 
haviors observed for different experimental systems. 

An interesting challenge will be to formulate experimental systems for which the 
ratio 8/E, can be finely tuned to check this hypothesis. 

• Most of the experimental results and theoretical analysis have dealt situations 
where the lifetime of bonds is large with respect to the relaxation time of the poly- 
meric linkers and so, dominates the systems response. This corresponds to brittle frac- 
tures: the understanding of the physics involved in nucleation and of brittle fractures 
begin to be clarified. The opposite situation where the relaxation of the polymeric 
linkers is comparable or larger than lifetime of the bonds is still poorly understood. 
Transient networks with polymer links in the entangled regime, or with nodes of tun- 
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able morphology, or near the percolation transition (marginal transient networks) will 
provide experimental realizations of this situation. From the theoretical side, this new 
time scale could be included in the SH-ABR model. This future direction will tackle 
dissipation mechanisms (viscous dissipation, structural reorganization, etc.) involved 
during nucleation of ductile fractures and during fracture propagation, to progress in 
the fundamental understanding of the concepts of brittleness and ductility in complex 
fluid. 

• Finally, the role of the rigidity of the polymeric linkers in the fracture of tran- 
sient networks is a big and challenging question to address in the future. 
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